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Introduction 



■ Recall that the reduction of a finite list of k, say, distinct (real) numbers [a±, 0,2, ■ ■ ■ , a^] is the 

unique permutation a = [a\, . . . , er^], of {1, . . . , k} such that a± is the ai-th largest element in the 
list, 02 is the <72-th largest element in the list, etc. In other words [0,1,0,2,... , a^] and a are "order- 
isomorphic". For example, the reduction of [6, 3, 8, 2] is [3, 2, 4, 1] and the reduction of [ir, 7, e, 0] is 
[4,1,3,2]. 



Given a permutation tt = ir\ ...7r n and another permutation c = [<ti, . . . ,0"^] (called a pattern), 
we denote by N a (Tr) the number of instances 1 < i\ < . . . < i& < n such that the reduction of 
7^ . . . vr ifc is a. 

For example, if tt = 51324 then 

-^[1,2,3] (tt) = 2 (because 7T27r37r5 = 134 and ^2^4,715 = 124 reduce to [1,2,3]). 
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■^[1,3,2] (tt) = 1 (because 7^37^ = 132 reduces to [1,3,2]). 
^[2,1,3] ( 7r ) = 1 (because ^37^5 = 324 reduces to [2, 1,3]). 

^[2,3,1] ( 7r ) = (because none of the 10 length-three subsequences of n reduces to 231). 

iV[ 3)1) 2](7r) = 5 (because 7ri7T27r 3 = 513 and 7ri7r 2 7T4 = 512 and 7ri7r 2 7r5 = 514 and 7ri7r 3 7T5 = 534 and 
7ri7T47r5 = 524 all reduce to [3, 1, 2]) . 

-W[3,2,i] (tt) = 1 (because 7ri7r 3 7T4 = 532 reduces to [3,2,1]). 

Of course the sum of iV CT (7r) over all /c-permutations a is (^). 

Fixing a pattern a, the set of permutations ir for which N^^n) = (we say that tt avoids a) is called 
the Wilf class of a, and more generally, given a set of patterns S, the set of permutations for which 
N cr (7r) = for all a € S, is the Wilf class of that set. The first systematic study of enumerating 
Wilf classes was undertaken in the pioneering paper by Rodica Simion and Frank Schmidt [SiSc]. 

The general question is extremely difficult (see [Wiki] and [Bo3]) and "explicit" answers are only 
known for few short patterns (and sets of patterns), the increasing patterns [1, 2, . . . , k], and a few 
other West- equivalent to them, giving the same enumeration. For example, even for the pattern 
[1, 3, 2, 4] (http : //oeis . org/A061552) the best known algorithm takes exponential time in n, and 
it is very possible that that's the best that one can do. 

But for those patterns a for which we know how to enumerate their Wilf classes, most importantly 
the increasing patterns [1, . . . , k], it makes sense to ask the more general question: 

Given a pattern a, and a positive integer r, find a "formula", or at least a polynomial-time algorithm 
(thus answering the question in the sense of Herb Wilf[Wil]) that inputs a positive integer n and 
outputs the number of permutations tt of {1, . . . , n} for which N a (7r) = r. We call such a class a 
generalized Wilf class. 

Ideally, we would like to have, given a pattern a, an explicit formula, in n and q, for the generating 
function (S n denotes the set of permutations of {1, . . . , n}) 

A a (q,n):= £ q N ^ j 

then, for any fixed r, the sequence of coefficients of q r in A a (q, n) would give the sequence enumer- 
ating permutations with exactly r occurrences of the pattern a. 

In fact, for patterns of length < 2 there are nice answers. Trivially 

A {1] {q,n) := n\q n , 
and almost-trivially (or at least classically) 

A [2jl] (g,n):=(l)(l + ( ? )...(l + g + ... + ^- 1 ) = [n}\ , 
the famous "(/-analog" of n\. But things start to get complicated for patterns of length 3. 
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Past Work 



For a very lucid and extremely engaging introduction to the subject, as well as the state-of-the-art 
in 2004, we strongly recommend Miklos Bona's masterpiece [Bo3]. 

In [NZ] , John Noonan and the second-named author initiated a functional equations-based approach 
for enumerating generalized Wilf classes. In order to illustrate it, they reproved John Noonan's[N] 
combinatorially-proved result that the number of permutations of length n with exactly one occur- 
rence of the pattern [1, 2, 3] equals ^ ( n ™ 3 ) • Recently a proof from the book of this result was given 
by Alexander Burstein[Bu] (see also [Zl]). 

In [NZ] it was conjectured that the number of permutations of length n with two occurrences of 
the pattern [1,2,3] equals |^2n-i)"n+5)° (n-J • This conjecture was proved by Markus FulmekfF], 
using Dyck paths. 

In [NZ] it was also conjectured that the number of permutations of length n with one occurrence 
of the pattern [1,3,2] equals n ~i)- This conjecture was proved by Miklos B6na[Bol], who 

later proceeded to prove[Bo2] the interesting fact that the sequences enumerating permutations 
with exactly r occurrences of [1,3,2] is P-recursive (i.e. satisfies a homogeneous linear recurrence 
with polynomial coefficients) for every r. In fact he proved the stronger result that the generating 
functions are always algebraic. This was vindicated by Toufik Mansour and Alek Vainshtein[MV] 
who gave an efficient algorithm to actually compute these generating functions, and they used it 
to find explicit expressions for 1 < r < 5. 

Another interesting but different "functional equation" approach, for patterns of length three, was 
developed by Firro and Mansour [FM]. 

This Project 

But so far, practically nothing is known for patterns of length larger than three and r > 0. In this 
paper we will modify the approach of [NZ] in order to generate, in polynomial time, such sequences 
for increasing patterns of any length [l,...,k]. That method can be extended to the patterns 
[1, . . . , k — 2, k, k — 1] and possibly other families, but here we will only discuss increasing patterns. 

Using the new algorithm to compute sufficiently many terms, we were able to conjecture explicit 
formulas, in n, for the number of permutations of length n with exactly r occurrences of the 
pattern [1,2,3], for 5 < r < 7, extending Fulmek's[F] conjectures for r = 3 and r = 4. We 
believe that the enumeration schemes, that our algorithms generate, should enable our computers 
to conjecture holonomic representations for the more general quantities (see below), that once 
guessed, should be amenable to automatic rigorous proving in the holonomic paradigm[Z2], using 
Christoph Koutschan's[K] far-reaching extensions and powerful implementations. But since these 
conjectures are certainly true, and their formal proof would (probably) not yield any new insight, 
we don't think that it is worth the trouble to actually carry out the gory details, wasting both 
humans' time (it would require quite a bit of daunting programming) and the computers' time (it 



3 



would take a very long time, due to the complexity of the schemes). 
Now let's recall the Noonan-Zeilberger Functional Equation Approach. 
The Noonan-Zeilberger Functional Equation Approach 

The starting point of the Noonan-Zeilberger [NZ] approach for enumerating generalized Wilf classes 
is to derive a functional equation. Let's review it with the simplest non-trivial case, that of the 
length-3 increasing pattern [1,2,3]. 

In addition to the variable q, introduce n extra catalytic variables x\, . . . , x n , and define the weight 
of a permutation ir = 7Ti . . . 7r n of length n by 

n 

weight^) := q N l^M x [^<-< b <^=^}\ , 

i=i 

(as usual, for any set A, \A\ denotes the number of elements of A). For example, 

weight(l23A5) = q w x\xlxlx± , 

weight{WZ21) = 1 , 
weight(21354) = q 4 x%xfxl = q 4 xfx^xl . 

Let's define the polynomial in the n + 1 variables 

P n (q;xi, . . . ,x n ) := ^ weight^) . 

Let 7r = 7Ti . . . 7r„ be a typical permutation of length n. Suppose 7Ti = i. Note that the number 
of occurrences of the pattern [1,2,3] in ir equals the number of occurrences of that pattern in the 
beheaded permutation 7r 2 . . . 7r n plus the number of the patterns [1, 2] in the beheaded permutation 
7T2 . . . 7r„ where the "1" is i + 1, or i + 2, or ... or n. Let ir' be the reduction to {1, . . . , n — 1} of 
that beheaded permutation. We see that 

weight(ir) = x™ % weight(ix ) \ Xi ^ qXi+l i X i +1 ^qx i+2 , ... , x n -i^>qx n 

The factor of x™~ 1 is because converting ir' from a permutation of {1, . . . , n — 1} to a permutation 
of {1, . . . , i — 1, i + 1, . . . , n}, and sticking an i at the front introduces n — i new [1, 2] patterns where 
the "1" is i. This gives the Noonan-Zeilberger Functional Equation 

n 

P n (q;xi, ...,x n ) = ^x"" l P n _i(^;xi, . . . , qx i+1 , . . .,qx n ) . (NZFEl) 

i=l 

Having found P n (q; xi, . . . , x n ), we set the "catalytic" variables xi, . . . , x n all to 1 and get 

f n {q) := A [h2 ,3](q,n) = P n (q; 1, 1, . . . , 1) . 
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Even though this is an "exponential-time" (and memory!) algorithm, it is much faster than the 
direct weighted counting of all the n! permutations, and we were able to explicitly compute them 
through n = 20. 

This is implemented in procedure fn(n,q) in P123. Procedure L20(q) ; gives the pre-computed sequence of fn(n,q) 
for n between 1 and 20 . 

Here are the first few terms: 

/i(?) = l , / 2 (?) = 2 , f 3 (q)=q + 5 , f 4 (q) = q 4 + 3 q 2 + 6 q + 14 , 



f 5 {q) =q W + 4q 7 + 6q 5 +9q 4 + 7q 3 + 24q 2 +27q + A2 , 



f 6 (q) = q 20 +5q 16 +8q 13 +6q 12 +6q 11 +16q w +12q 9 +24q 8 +32q 7 +37q 6 +54q 5 +74:q 4 +70q 3 +133q 2 +n0q+132 , 



f 7 (q) = q 35 + Q g 30 +10 q w + 10 q 25 +8 q 23 + 13 q 22 +30 q 21 + W q 20 +32 g 19 + 18 g 18 +62 q 17 +7A g 16 +24 g 15 + 100 q 14 
+130 g 13 +104 g 12 +162 g n +191 <? 10 +232 g 9 +260 <? 8 +320 g 7 +387 <? 6 +395 g 5 +507 <? 4 +461 <? 3 +635 g 2 +429 g+429 , 



/ 8 (q) =q 56 + 7 q 50 + 12 q 45 + 15 q 44 + 10 q 41 + 16 q 40 + 40 q 39 + 18 q 38 + 47 q 36 + 38 q 35 + 68 q 34 + 60 q 33 

+58 q 32 + 66 g 31 + 154 q 30 + 138 q 29 + 1 15 q 28 + 156 q 27 + 252 q 26 + 324 q 25 + 228 q 24 + 288 q 23 + 537 q 22 
+466 q 21 + 546 q 20 + 656 q 19 + 682 q 18 + 1004 q 17 + 1047 g 16 + 886 q 15 + 1494 g 14 + 1456 q 13 + 1580 q 12 
+1818 g n +2077 g 10 +2182 g 9 +2389 g 8 +2544 g 7 +2864 g 6 +2570 g 5 +3008 g 4 +2528 g 3 +2807 g 2 +1638 g+1430 

For f n (q) for 9 < n < 20 see: 

http : //www. math. rutgers . edu/~zeilberg/tokhniot/oP123d . 

Using this data, the computer easily finds rigorously-proved explicit expressions for the first six 
moments (about the mean) of the random variable "number of occurrences of the pattern [1, 2, 3]", 
and from them verifies that, at least up to the sixth moment, this random variable is asymptocally 
normal, as humanly proved (for all patterns) by Miklos B6na[Bo4]. See: 

http : //www. math. rutgers . edu/~zeilberg/tokhniot/oP123a . 
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The "Perturbation" Approach 



The equations of quantum field theory are (usually) impossible to solve exactly, but physicists got 
around it by devising clever "approximate" methods using perturbation expansions, that only use 
the first few terms in a potentially "infinite" (and intractable) series, but that suffice for all practical 
purposes, using Feynman diagrams. 

Of course, we are enumerators, and we want exact results, but suppose we only want to know the 
sequences enumerating permutations with exactly s occurrences of the pattern [1, 2, 3] for s < r for 
some relatively small r, rather than for r = (™), provided by the full 

f n (q) = P n (q; 1 [n times}) . 

In the original article [NZ], for r = 0, Noonan and Zeilberger simply plugged-in q = and x\ = 
. . . = x n = 1, getting a simple enumeration scheme, that proved, for the n-th time, the classical 
result that the number of permutations of length n that avoid the pattern 123 equals the Catalan 
number (2n)!/(n!(n + 1)!). For r = 1, they differentiated Eq. (NZFEl) with respect to q, using 
the multivariable calculus chain rule, and then plugged-in q = and x\ = . . . = x n = 1. For r = 2 
they did it again, but this turned out to be, for larger r, a Rube Goldberg nightmare, even for a 
computer. 

Here is a much easier way! 

Recall that you are really only interested in f n (q) = P n (q; 1 [n times] ). Plugging it into (NZFEl) 
gives 

n 

P n (q; 1 [n times ]) = P n _i{q ; 1 [i — 1 times ] , q[n — i times]) 

i=i 

This forces us to put-up with expressions of the form 

P ao+ai (q; 1 [ao times], q[a x times]) . 
Plugging this into (NZFEl) yields 

P ao+ai (<?; 1 [a times] , q[ai times]) = ^P ao+ai _i(g; 1 [i-l times] , q[a -i times] , q 2 [ai times]) 

ai 

+ g ai ~*P an+ai _i(g; 1 [a times] , q [ i — 1 times] , q 2 [ a\ — i times]) 

i=l 

This forces us, in turn, to consider expressions of the form 

P ao+ai+a2 (q; l[a times], q[a x times], q 2 [a 2 times]), 
that would force us to further consider expressions of the form 

P an+ai+a2+a3 (q; 1 [a times], q[a x times], q 2 [a 2 times] , q 3 [a 3 times]) , 
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etc. etc., leading to an exponential explosion in both time and memory. 

But, if we are only interested in the first r coefficients of f n (q), then we can take advantage of the 
trivial (you prove it!) but crucial lemma. 

Crucial Lemma: For s > r + f , the coefficients of q°, q 1 , . . . , q r of 

P ao+ai+ ... +as (q; l[a times ],... a s_1 ,[ a s _i times ], q s [ a s times]) 

- P ao+ai +...+a s (q; 1 [ a times] , . . . , q r [a r times ] , q r+1 [ a r +i + a r+2 + . . . + a s times]) 

all vanish. In other words, if for any polynomial p(q) in q, p^ (q) denotes the polynomial of degree r 
obtained by ignoring all powers of q larger than r, and letting CHOP r [p(q) ] := p^ r \q), (NZFE1) 
becomes (below, let n := a\ + . . . + a r + a r +i, and for any expression R and positive integer k, 
R$k, denotes R . . . R[k times] for example q 2 $3 means q 2 ,q 2 , q 2 ) 

P^\q; f$a , q%ai , ... , q r $a r ,q r+1 $a r+1 ) 

ao 

= CHOP r (^2pi r } 1 (q; - 1) , q$(a - i) , q 2 $ ai , g r $a r _! , q r+1 $(a r + a r+1 )) 

i=l 

ai 

+ q ai - l+a2+ - +a ^P^ r } 1 (q; l$a , q$(i - 1) , g 2 $(ax - i) , . . . , ^$a r _x , q r+1 $(a r + a r+1 )) 

i=i 

+ £g2(a a -i+as+...+a r+1 )pW i((z . ^ ? ^ ? g 2 $( g 3 $(Q2 _ •) > _ _ ^ (f $a r _! , <f + 1 $(a r .+a r .+l)) 
i=l 

+ 

ar+i 

+ q (r+D("r + i-i)p(r) i{q . 1$flo ? ? q 2 $(l2 ^ q r Ur ^+1$^ _ 1))) . 

i=l 

Now note that, because of the CHOP r operator in front, many terms automatically disappear, 
because of the powers of q in front. The bottom line is that the computer can automatically 
generate a scheme for computing the degree-r polynomials in q, 

F r (a ,...,a r+1 )(q) := Pj£+... +ar+1 (q; 1 [ao times], q [a,! times] q r+1 [a r+1 times]) , 

with ao + . . . + a r +i = n and ao, • • • , a r +i > 0. The number of such quantities is the coefficient of z n 
in 1/(1 - z) r+2 that equals (-l) r+2 (~ (r n +2) ) = { r+ r +V) terms - So each iteration involves 0(n r+1 ) 
evaluations and hence 0(n r+2 ) additions and doing it n times yields an 0(n r+3 ) algorithm for 
finding our object of desire, the degree r polynomial in q: 

fP(q) = F r (n,0[r + 1 times]){q) . 

Having found the scheme, the very same computer (or a different one), may use it to generate as 
many terms as desired. 
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The Maple package PI 23 



The Maple package P123 downloadable from 

http : //www. math. rutgers . edu/~zeilberg/tokhniot/P123 , 

implements the functional equation (NZFE1) and easily generated the first 25 terms of the enu- 
merating sequences for < r < 7. With this data, it empirically verified the already-known results 
for the number of permutations with exactly r occurrences of the pattern [1,2,3] for < r < 2, 
and made conjectures for r < 7 as follows. Let a r {n) be the number of permutations of length n 
with exactly r occurrences of the pattern [1,2,3]. 

/ \ n (2n-l)! 
a (n) = 2 



ai(rt) = 6 



(n - 1)! (n + 1)! 

(2n - 1)! 
(n-3)! (n + 3)! 



(2n — 2V 

^")= (n -4) !( n + 5)r < 59 " 2 + 117 " + 100) 



a 3 (n) = t — ^ 2n , . 3 ^ ! w ., • An (113 n 3 + 506 n 2 + 937 n + 1804) 
(n — 5)! (n + 7)! ' 



(2n-4)! 

a4(n) - 



(n-4)!(n + 9)! 

(3561 n 8 + 3126 n 7 - 46806 n 6 + 12384 n 5 - 659091 n 4 + 2630634 n 3 + 5520576 n 2 + 26283456 n - 39191040) 



(2n-5)! 

a 5 (n) = 



(n-5)!(n + ll)! 

( 26246 n 10 + 136646 n 9 - 115872 n 8 + 22524 n 7 - 9648450 n 6 + 71304534 n 5 
+381205612 n 4 + 1607633896 n 3 + 2800103664 n 2 + 3611692800 n - 32891443200) 



(2n-6)! 

a 6 (n) - 



(n-6)!(n + 13)! 

( 193311 n 12 + 2349954 n 11 + 13035003 n 10 + 95151030 n 9 + 406430793 n 8 + 2889552582 n 7 
+14335663329 n 6 + 60005854890 n 5 + 313010684796 n 4 + 1025692693464 n 3 
+1283595375168 n 2 - 6909513045120 n - 28177269120000 ) . 
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(2n - 7)! 
° 7 W- ( n -5)!(n + 15)!' 

( 1386032 n 16 + 13111080 ra 15 + 22526480 n 14 + 355187760 n 13 - 1654450096 n 12 + 10534951680 n 11 
+15797223760 n w - 305671694640 n 9 + 3750695521216 n 8 - 26631101348520 n 7 
-86395090065440 n 6 - 636425872408320 n 5 + 3647384624274048 n 4 
+11386434230674560 n 3 +103032675524966400 n 2 -157858417817856000 n-763734137886720000 ) 

Enumerating Permutations with r occurrences of the pattern [1,2,3,4] for small r via 
a Noonan-Zeilberger Functional Equation 

In addition to the variable q, we now introduce 2n extra catalytic variables x±, . . . , x n , and j/i, . . . , y n , 
and define the weight of a permutation ir = 7Ti . . . 7r n of length n by 

n 

Weight^) := ^[1,2,3,4] M FJ a ,l{ 1 < a < 6 < c <™ ; 7r a =i<7n,<7r c }| _ y\{l<a<b<n ; w a =i<-K b }\ 

i=l 

For example, 

weight([l, 2, 3, 4, 5, 6]) = g 1 ^ ^^!^!^ , 
weight([6,5,4,3,2, 1]) = 1 , 
weight([3, 4, 5, 6, 1]) = qx\x^y\y\y^ . 



Let's define the polynomial in the 2n + 1 variables 

P n (^;xi,...,x n ;yi,...,y n ) := \ ] weight^) 



E 



Let 7r = 7Ti . . . , 7r n be a typical permutation of length n. Suppose 7Ti = i. Note that the number 
of occurrences of the pattern [1,2,3,4] in ir equals the number of occurrences of that pattern 
in the beheaded permutation ir 2 ■ ■ ■ ir n plus the number of the patterns [1,2,3] in the beheaded 
permutation 7r 2 . . . 7r n where the "1" is i + 1, or i + 2, or . . . or n. Let it' be the reduction to 
{1, . . . ,n — 1} of that beheaded permutation. Also note that the number of occurrences of the 
pattern [1,2,3] where the "1" is an i gets increased by the number of occurrences of the pattern 
[1, 2] in the beheaded permutation, where the "1" is a j with j > i. We see that 

weight{ir) = y i weight(ir ) \ Xi ^ qXi+l i Xi+l ^ qXi+2t ... i Xn _ 1 ^ qXn ■ yi ^ Xiyi+1 i y i+1 ^xiy i+2 , ... , y n -i^xiy n 

The factor of is because converting it' from a permutation of {1, . . . , n — 1} to a permutation 
of {1, . . . , i — l,i + l,...,n}, and sticking an i at the front introduces n — i new [1, 2] patterns where 
the "1" is i. This gives the Noonan-Zeilberger Functional Equation for the pattern [1,2,3,4]: 

n 

P n (q;xi,. ..,x n ;y 1 ,...,y n ) = ^ y"" l P„_i(g; x u . . . , Xi-i, qx i+1 , . . . , qx n ; y u . . . ,yi-i,Xiy i+ i, . . .,Xiy n ) 



i=l 



(NZFE2) 
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Having found P n (q; xi, . . . , x n ; yi, . . . , y n ), we set the "catalytic" variables x\,...,x n and yi, ■ ■ ■ ,y n 
all to 1 and get 

g n {q) := A [lj2 ,3,4](q,n) = P n (q; 1, 1, . . . , 1 ; 1,1,... ,1) . 

Even though this is an "exponential-time" (and memory!) algorithm, it is still faster than the 
direct weighted counting of all the n! permutations, and we were able to explicitly compute them 
through n = 10. 

The first few polynomials are 

9i(q) = 1 , 92{q) = 2 , g 3 (q) = 3 , g 4 (q) = q + 23 , 
9 5 (q) =q b + Aq 2 + 12^ + 103 , 

5e(g) =q 15 + 5q 9 + 8q 6 + 12q 5 + 6 g 4 + 10 g 3 + 63 g 2 + 102^ + 513 , 



g 7 (q) = q 35 + 6q 25 + 10 q 19 + 18 q w + 12 g 15 + 13 g 13 + 24q u + 32q 10 +72q 9 + 10 g 8 + 46<? 7 
+142 g 6 + 116 q 5 + 146 g 4 + 196 q 3 + 665 g 2 + 770 q + 2761 , 



58 ( g ) = g 70 +7 g 55 + 12 g 45 + 15 g 41 + 10 q 39 +8 q 36 +28 q 35 +40 q 32 +41 q 29 + W q 2S +24 q 27 +U q 26 +84 q 25 

+24 q 2A + 89 q 23 + 12 q 21 + 142 q 20 + 136 q 19 + 96 g 18 + 115 q 17 + 333 q 16 + 156 q 15 + 112 g 14 + 312 q 13 

+199 g 12 +600 <? n +573 g 10 +804 g 9 +503 g 8 +885 <? 7 +1782 <? 6 +1204 g 5 +2148 <? 4 +2477 g 3 +5982 g 2 +5545 g+15767 

For gg (q), gio(q) see: http://www.math.rutgers.edu/~zeilberg/tokhniot/oP1234d. 

The obvious analog of the Crucial Lemma still holds, and one can get polynomial time (in n) 
algorithms, to compute the number of permutations of length n with exactly r occurrences of 
the the pattern [1,2,3,4]. Alas, because we have twice as many catalytic variables, the 0{n r+3 ) 
becomes 0{n 2r+b ). Nevertheless, we were able to compute the first 70 terms for the case r = 1. 
Here are the first 23 terms: 



0, 0, 0, 1, 12, 102, 770, 5545, 39220, 276144, 1948212, 13817680, 98679990, 
710108396, 5150076076, 37641647410, 277202062666, 2056218941678, 15358296210724, 
115469557503753, 873561194459596, 6647760790457218, 50871527629923754 . 
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The rest can be viewed in: http://www.math.rutgers.edu/~zeilberg/tokhniot/oF1234a. 
The Maple package PI 234 

Everything is implemented in the Maple package PI 234. See the webpage of this article for sample 
input and output files. The package F1234 is a more efficient implementation for a small number 
of occurrences r. 

Beyond 

Of course, the same reasoning applies to any increasing pattern [1, . . . , k] but we have, in addition 
to q, (k — 2)n additional catalytic variables. For each specific r this implies a scheme that enables 
one to compute in "polynomial" time (in n, but of course not in k or r) the desired numbers. For 
the patterns [1,2,3,4,5] and [1,2,3,4,5,6] (i.e. k = 5 and k = 6) this is implemented in Maple 
packages P12345 (and its more efficient [for small r] version F12345) and P123456 respectively. 

Other Patterns 

Even the case of pattern-avoidance, i.e. r = 0, is already extremely difficult in general. As we 
mentioned above for the pattern [1,3,2,4], there is no known polynomial time algorithm for enu- 
merating permutations that avoid it. But for some few infinite families (see [Bo3] and [Wiki]) exact 
formulas (for the avoiding, r = 0, case) are known, and the present approach would hopefully be 
able to find polynomial-time schemes for r > 0, at least for some of them. We hope to investigate 
this in a future paper. 
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